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Abstract 

Using recent precise hadronic r-decay data on the V — A spec- 
tral function, and general properties of QCD such as analyticity, the 
operator product expansion and chiral perturbation theory, we get 
accurate values for the QCD chiral order parameters L\q{Mp) and 
Cgj{Mp). These two low-energy constants appear at order and p^, 
respectively, in the chiral perturbation theory expansion of the V — A 
correlator. At order we obtain L^o(^p) = -(5-22 ± 0.06) • lO'^. 
Including in the analysis the two-loop (order p^) contributions, we 
get L[q{Mp) = -(4.06 ± 0.39) • lO^^ and Cl^{Mp) = (4.89 ± 0.19) • 
10~^ GeV~^. In the SU(2) chiral effective theory, the corresponding 
low-energy coupling takes the value Is = 13.30 it 0.11 at order p^, and 
I5 = 12.24 ± 0.21 at order /. 



1 Introduction 



The precise hadronic r-decay data provided in refs. [1-6] are a very impor- 
tant source of information, both on perturbative and non-perturbative QCD 
parameters. The theoretical analysis of the inclusive r decay width into 
hadrons allows to perform an accurate determination of the QCD coupling 
as{Mr) [7-11], which becomes the most precise determination of as{Mz) af- 
ter QCD running. In this case, non-pcrturbative QCD effects parametrised 
by power corrections are strongly suppressed. Another example of the use 
of hadronic r-decay data is the study of SU(3)-breaking corrections to the 
strangeness-changing two-point functions [12-16]. The separate measure- 
ment of the IAS*! = and \AS\ = 1 tau decay widths provides accurate 
determinations of fundamental parameters of the Standard Model, such as 
the strange quark mass and the Cabibbo-Kobayashi-Maskawa quark-mixing 
IK.I [16]. 

Very important phenomenological hadronic matrix elements and non- 
perturbativc QCD qiiantities can also be obtained from r-dccay data. Of 
special interest is the difference of the vector and axial-vector spectral func- 
tions, because in the chiral limit the corresponding V—A correlator is exactly 
zero in perturbation theory. The r-decay measurement of the V — A spec- 
tral function has been used to perform [17-19] phenomenological tests of the 
so-called Weinberg sum rules (WSRs) [20], to compute the electromagnetic 
mass difference between the charged and neutral pions [18], and to determine 
several QCD vacuum condensates [21,22]. Prom the same spectral function 
one can also determine the A7 — 3/2 contribution of the A^" = 1 four-quark 
operators and Qs to s'^/ek-, in the chiral limit [23]. 

Using chiral perturbation theory (xPT) [24-26], the hadronic r-decay 
data can also be related to order parameters of the spontaneous chiral sym- 
metry breaking (SxSB) of QCD [27]. xPT is the effective field theory of QCD 
at very low energies; it describes the SxB Nambu-Goldstone boson physics 
through an expansion in external momenta and quark masses. The coeffi- 
cients of that expansion are related to order parameters of S^SB. At lowest 
order (LO), i.e. O(p^), all low-energy observables are described in terms of 
the pion decay constant ~ 92.4 MeV and the light quark condensate. 
At next-to-leading order (NLO), 0{p^), the SU(3) xPT Lagrangian contains 
12 low-energy constants (LECs), Lj=i,...,io and Hi^2 [26]. At 90 (23) 

additional parameters Ci=i,... ,90 appear in the even (odd) intrinsic parity sec- 
tor [28] . These LECs are not fixed by symmetry requirements alone and have 
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to be determined phenomenologically or using non-perturbative techniques. 
The 0{p'^) Li couphngs have been determined in the past to an acceptable 
accuracy; a recent compilation can be found in ref. [29]. Much less well 
determined are the C(p^) couplings Q. 

There has been a lot of recent activity to determine the chiral LECs from 
theory, using as much as possible QCD information [30-39]. This strong 
effort is motivated by the precision required in present phenomenological ap- 
plications, which makes necessary to include corrections of 0{p^). The huge 
number of unknown couplings is the major source of theoretical uncertainty. 

In this paper we present an accurate determination of the xPT couplings 
Lio and Cs7, using the most recent experimental data on hadronic r decays 
[1]. Previous work on Liq using r-decay data can be found in refs. [18,19, 
21,40]. Our analysis is the first one which includes the known two- loop xPT 
contributions and, therefore, provides also the 0{p^) coupling Csr- 



2 Theoretical Framework 

The basic objects of the theoretical analysis are the two-point correlation 
functions of the vector and axial- vector quark currents, defined as follows: 

n^^(5) = ijd'xe^'^^o\T{ji;{x)jr,{oy)\o) 

(1) 

Here, we just need the non-strange correlators, i.e. J'ij{x) denotes the 
Cabibbo-allowed vector or axial- vector currents, V^^{x) = uj^d and A'^^ = 
uj'^^^d. Moreover, our analysis will concentrate in the difference 



+ n(s), (2) 



_ 2/1 



s — ml 



where we have made explicit the contribution of the pion pole to the lon- 
gitudinal axial-vector two-point function. We will work in the isospin limit 



TTiu = md where ^ulviq"^) = 0- 
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The correlator n(s) is analytic in the entire complex s-plane, except for 
a cut on the positive real axis which starts at the threshold Sth = ^m^. Ap- 
plying Cauchy's theorem to the circuit in Fig. [H one gets the exact relation: 

so I I r 
dss"-Imn(s) + (b dss"n(s) 

= 2f^ml- + Res[s'^U{s),s = 0]. (3) 

For non-negative values of the integer power n, the pion pole is the only 
singularity within the contour and one gets the so-called finite energy sum 
rules (FESR), widely used in the literature. When n takes negative values, 
the weight factor introduces a pole at the origin which gives rise to the 
additional contribution in the r.h.s. of the equation, given by the residue of 
s"n(s) at s = 0. 

In the chiral limit (m„ = = 0) the correlator n(s) vanishes identically 
to all orders in perturbation theory. For large enough Euclidean values of 
s = —Q^ its operator product expansion (OPE), n(Q^) = Yl,k^Xk^ IQ'^^ i 
contains only power-suppressed contributions from dimension d = 2k opera- 
tors, starting at = 6. The nonzero up and down quark masses induce tiny 
corrections with dimensions two and four, which are negligible at high values 
of Q^. Therefore, with n > and Sq large enough so that the OPE can be 
applied in the entire circle s = Sq, the integral over the spectral function from 
Sth to So is equal to the pion pole term plus the OPE contribution 

{—i)"'CY,~^-.s generated by the integration along the circle. For n = and 



2(n+l) 

n = 1, C*^~^-^-| is zero in the chiral limit and one gets the celebrated first and 
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second WSRs [20], respectively. 

For negative values of n = —m < 0, the OPE does not give any contri- 
bution to the integration along the circle s = sq. One gets then: 



ds 1 



TT 



Imn(s) 



2m 
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(m — 1) 



n(™-i)(o) 



(4) 



(m- 1)! 

where 11 (0) denotes the (m — l)th derivative of n(s) at s = 0. The 
interest of this relation stems from the fact that at low values of s the cor- 
relator can be rigourously calculated within xPT. At present n(s) is known 
to 0{p^) [41], in terms of the LECs that we want to determine. The choices 
m — 1 and m — 2 allow then us to relate the spectral function measured in 
T decays with the theoretical expressions of n(0) and n'(0), which can be 
derived from the results obtained in ref. [41]: 
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where Hi = mj\og{mi/ n) / {I&t:'^ fl). 

To a first approximation the effective parameters L^q and correspond 
to the LECs L[q{ii) and C^'^{ii), respectively. At O(p^), the only relevant 



4 



correction is given by the logarithmic terms in the second hne of which 
cancel the xPT renormalization scale dependence of L\q{^); these contribu- 
tions are suppressed by one power of l/Nc with respect to L^Q^fi), where Nc 
is the number of quark colours. The rest of lines in ([5]) contain the 0{p^) 
corrections: the tree-level contributions from the 0{p^) xPT Lagrangian are 
given in the third and fourth lines, the term proportional to (Lg + 2L^q)(//) 
in the fifth line is the one- loop contribution of the 0{p'^) xPT Lagrangian, 
and the function G2L(At, s = 0) in the last line, which does not depend on any 
LEG, contains the proper two-loop contributions. 

In Eq. ([6]) the tree-level contribution is given by Cgj{fi), whereas the 
term proportional to L^fi) is a one- loop correction, which is suppressed 
by one power of I/Nq, and the two-loop contributions are contained in 
^2l(/^5 s) = ^G'2l(/^5 s). The derivative operation, when acting over the one- 
loop contribution to n(s), generates the terms proportional to inverse powers 
of the pion and kaon masses in the second line. For simplicity, we omit the 
explicit analytic forms of G2L{fJ') and G'2^(/i), which are very lengthy and 
not too enlightening; these two functions contain a 1/Nq suppression factor 
with respect to L^Q^fi) and C^^^fi). 

3 Determination of Effective Couplings 

We will use the 2005 ALEPH data on semileptonic r decays [1], which pro- 
vides the most recent and precise measurement of the V —A spectral function. 
The effective chiral couplings can be directly extracted from the following in- 
tegrals over the hadronic spectrum: 

-8L5|^ ^ n(0) = - / -Imn(.), (7) 

TT /„ , S 

leCgt = n'(o) = - / -imn(s). (8) 

These relations are exactly satisfied at sq ^ oo. At finite values of sq, they 
assume that the OPE approximates well the correlator n(s) over the entire 
complex circle |s| = sq. The OPE is expected to be a valid approximation 

^Or equivalently these relations assume that the integrals on the real axis from sq to 
infinite are negligible, what is expected to be true only for high enough values of sq and 
for accidental "duality points" . 
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for high-enough values of Sq and away from the real axis. While the kinemat- 
ics of T decay restrict the upper limit of integration to the range Sq < m^, the 
main source of theoretical uncertainty in the contour integration originates in 
the region close to the point s = sq in the real axis. Studying the sensitivity 
to So of the integrals ([7j) and ([8]), one can test validity of the OPE and assess 
the size of the associated systematic errors. 




-0.01 ot ' / 



Figure 2: Determinations of LIq at different values of sq. The continuous 
lines show the results obtained from Eq. ([7]). The modified expressions in 
Eqs. ([9]) and (fTOj) give rise to the dashed and dot-dashed lines, respectively. 
For clarity, we do not include their corresponding error bands. 

In Fig. [21 we plot the value of L^q obtained from Eq. ([7]) for different 
values of sq. The band between the continuous lines shows the corresponding 
experimental uncertainties (at one sigma). As expected, the result is far from 
an horizontal line at low values of sq, where the applicability of the OPE 
is suspect. The oscillatory behaviour stabilises quite fast reaching a rather 
stable and flat result at values of Sq between 2 and 3 GeV^. The weight factor 
1/s decreases the impact of the high-energy region, minimising the size of 
quark-hadron duality violations around sq. This integral appears then to be 
much better behaved than the corresponding FESRs with s" (n > 0) weights. 

There are several possible strategies to estimate the central value for L^q 
and the unavoidable theoretical uncertainties. One is to give the predictions 
fixing So at the so-called "duality points" , where the first and second WSRs 
happen to be satisfied. Owing to the oscillatory behaviour of the WSRs 
results, this happens at two different values of sq. At the highest "duality 
point" , which is obviously the more rehable, we obtain Ll^ = -(6.45 ±0.09) • 
10"'^, where the quoted error only includes the experimental uncertainty. 
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Being very conservative, one could also take into account the first "duality 
point"; performing a weighted average of both results, we get L^q = — (6.50± 
0.13) ■ 10"'^, where the uncertainty covers the values obtained at the two 
"duality points" . 

Assuming that the integral ([7]) oscillates around his asymptotic value 
with decreasing oscillations, one can get another estimate performing an 
average between the maxima and minima of the successive oscillations. This 
procedure gives a value L^g = —(6.5 ±0.2) ■ 10~^, that is perfectly compatible 
with the previous results based on the "duality points". Our last method 
of estimating the quark-hadron duality violation uses appropriate oscillating 
functions defined in [42] which mimic the real quark-hadron oscillations above 
the data. These functions are defined such that they match the data at 
approximately 3 GeV^, go to zero with decreasing oscillations and satisfy the 
first and second WSRs. We find in this way Lfl = -(6.50 ± 0.12) ■ lO'^, 
where the error spans the range generated by the different functions used. 
This result agrees well with our previous estimates. 

We can take advantage of the WSRs to construct modified sum rules with 
weight factors proportional to (1 — s/sq), in order to suppress numerically 
the role of the suspect region around s ~ Sq [8]: 



'Lfo = - r - fl--) Imn(s) + Ai(so), 



(9) 



so 



ds ( ^ s 



2 



— 1 Imn(s) 

rrJs,, s \ So J 

+ 2Ai(so) - A2(so). (10) 



The factors A,{so) = (2/^ + C^^^) / so and A^iso) = (2/>^ - Cj-^) /si 
are small corrections dominated by the term, since C2 4^ vanish in the 
chiral limit. The sum rule (fTOj) has been previously used in refs. [21,40]. 

The dashed and dot-dashed lines in Fig. [2] show the results obtained from 
Eqs. ([9]) and (ITOl) . respectively. As already found in refs. [21,40], the modified 
weight factors minimise the theoretical uncertainties in a very sizeable way, 
giving rise to very stable results over a quite wide range of sq values. One gets 
thenLl^ = -(6.51±0.06) • lO^^ using Eq. ©, and Ll^ = -(6.45±0.06) ■ 10~^ 
from Eq. (fTOl). 

Taking into account all the previous discussion, we quote as our final 
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result: 



-(6.48 ±0.06) ■ 10"^ 



(11) 
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Figure 3: Determinations of Cgy at different values of sq. The continuous 
lines show the results obtained from Eq. ([H]). The modified expressions in 
Eqs. (fT2|) and (fT3|) give rise to the dashed and dot-dashed lines, respectively. 
For clarity, we do not include their corresponding error bands. 

We have made a completely analogous analysis to determine the effective 
coupling . The results are shown in Fig. [31 The continuous lines, obtained 
from Eq. ([8]), are much more stable than the corresponding results for L'^, 
owing to the factor in the integrand. The discontinuous and dotted lines 
correspond to the results obtained from the modified sum rules: 
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The agreement among the different estimates is quite remarkable. We quote 
as our final conservative result, 
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•-^87 



(8.18 ±0.14) ■ lO^^GeV" 



(14) 
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4 Determination of L^q and Cg^ 

The xPT coupling L\q{h) can be obtained from L^q, using the relation ([5]). At 
0{p'^) the determination is straightforward, since one only needs to subtract 
from L^o the term [l — log (/i^/m^) + | log (m|-/m^)] /(1287r^). Taking /i = 
Mp as the reference value for the xPT renormalization scale, one gets 

L^q(Mp) = -(5.22 ±0.06) ■ 10~^ (15) 

At order p^, the numerical relation is more subtle because it gets small 
corrections from other LECs. It is useful to classify the 0{p^) contribu- 
tions through their ordering within the 1/Nc expansion. The tree- level term 
A.m\{Cl^ — ~ Cm){^p)^ which is the only 0{p^) correction in the large- 
Nc limit, is numerically small because it appears suppressed by a factor 
m^. The three relevant couplings have been determined phenomenologically 
with a moderate accuracy: Qi(Mp) = (1.24 ± 0.44) ■ 10"^ GeV"^ [34] (from 
ni'S;^(0) - ni°+^)(0)), CUM,) = (0.4 ± 6.3) ■ 10-5 GeV-2 [43] (from the 
Ktt scalar form factor) and C^soi^p) = (2-1 ± 0.5) ■ lO'^ GeV"^ [44] (from 
ai/Ki mass and width differences). These determinations agree reasonably 
well with published meson-exchange estimates [35, 41] and lead to a total 
contribution 4m2(Q^ - C^^ - Qo)(^p) = -(6.7 ± 5.2) • 10"^ The scale de- 
pendence of this combination of 0{p^) couplings [28] between fi = 0.6 GeV 
and fi = 1.1 GeV is within its quoted uncertainty. 

At NLO in 1/Nc we need to consider the tree-level contribution pro- 
portional to the combination of LECs (Cg2 — Cy^ — Cgi){Mp). We are not 
aware of any published estimate of these 1/Nc suppressed couplings, be- 
yond the trivial statement that they don't get any tree-level contribution 
from resonance exchange [35]. We will adopt the conservative range |Cg2 — 
- C'siK^p) < IQi - C'u - QoK^pi/S, which gives a contribution 
4(2m^ + m2)(Q2-C!3-C8i) W = (0.0±5.8) ■ 10^^ The scale dependence 
between /i = 0.6 GeV and = 1.1 GeV of this combination of 0{p^) cou- 
plings [28] is within its quoted uncertainty. The uncertainty on this term will 
dominate our final error on the L\q{Mp) determination. At the same NLO 
in 1/Nc, there is also a one- loop correction proportional to Ll{Mp); using 
the C(/) determination L^(Mp) = (5.93 ±0.43) ■ lO'^ [45], this contribution 
can be estimated to be 2(2/i^±^j^)L^(Mp) = -(1.56 ± 0.11) ■ 10"^ Finally, 
the 1/N^ suppressed two-loop function which collects the non-analytic con- 
tributions takes the value G2L{Mp) = —0.524 ■ 10"^, one order of magnitude 
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smaller than L^q, but still eight times larger than the uncertainty quoted for 
L^Q in fill I) . Taking all these contributions into account, we finally get the 
wanted 0{p^) result: 

L^o(Mp) = -(4.06±0.04ieff±0.39LECs)-10-=' 

= -(4.06 ± 0.39) ■ 10-^ (16) 

where the uncertainty has been split into its two main components. The 
final error is completely dominated by our ignorance on the 1 / Nc suppressed 
LECs of 0(/). 

The determination of Cg^ from Cgf does not involve any unknown LEG. 
The relation (jS]) contains a one-loop correction of size —(3.15 ± 0.13) ■ 10~^, 
which only depends on Lg(Mp) and the pion and kaon masses, and small non- 
analytic two-loop contributions collected in the term G'2i{Mp) = —0.277 • 
10~^ GeV~^. In spite of its 1/Nc suppression, the one-loop correction is very 
sizeable, decreasing the final value of the 0{p^) LEG: 

C;^{Mp) = (4.89 ± 0.19) ■ 10"^ GeV"2 . (17) 

5 SU(2)xPT 

Up to now, we have discussed the LEGs of the usual SU(3) xPT. It turns use- 
ful to consider also the effective low-energy theory with only two flavours of 
light quarks. In some cases, this allows to perform high-accuracy phenomeno- 
logical determinations of the corresponding LEGs at NLO. Moreover, recent 
lattice calculations with two dynamical quarks are already able to obtain the 
SU(2) LEGs with sufficient accuracy and this is an important check for them. 

In SU(2) xPT, there are ten LEGs, /i=i,..7 and /ii,2,3, at O(p^) (NLO) 
[25]. Using the 0{p^) relation between IKfi) and L5q(/x), recently obtained 
in ref. [46], and the definition of the invariant couplings k adopted in [25], 



10 



we get 
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(18) 



where rh\ = ml — ml/2 is the kaon mass squared in the limit m„ = = 0, 
Xk = m|-/(167rV^), f^K = rh\ log(mi^//i)/(167rV^) and pi ~ 1.41602. 

The first line contains the 0{p'^) contributions; the determination of 
at this order is then straightforward. The full C(p^) result, with the differ- 
ent tree-level, one-loop and two-loop corrections, is given in the other lines. 
Following the same procedure as in the SU(3) case, we get the results 



13.30 ±0.11, 
12.24 ±0.21, 



(19) 



6 Summary 

Using the most recent hadronic r-decay data [1] on the V — A spectral func- 
tion, and general properties of QCD such as analyticity, the OPE and xPT, 
we have determined very accurately the chiral LECs L'[q{Mp) and Clj{Mp). 
Performing an 0{p'^) analysis, we obtain 

Llo(Mp) = -(5.22 ± 0.06) • 10"^ , (20) 

while a more elaborate study, including the 0{p^) xPT corrections provides 
the values: 

L^o(Mp) = -(4.06±0.04icft ±0.39lecs)-10-' 

= -(4.06 ± 0.39) • 10-^ (21) 
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and 

C^^(Mp) = (4.89 ± 0.19) ■ 10-3 GeV"2 _ (22) 

Our error estimate includes a careful analysis of the theoretical uncertain- 
ties associated with the use of the OPE in the dangerous region close to the 
physical cut. Moreover, in (12T|) we have explicitly separated the error into its 
two main components, showing that our present ignorance on the I/Nq sup- 
pressed LECs dominates the final uncertainty of the L^q(Mp) determination 
at C(/). 

Several determinations of Lio have been performed before [18,19,40], 
using the older 1998 ALEPH data [2,3]. In ref [18] the result L^o(Mp) = 
— (5.13 ± 0.19) ■ 10~3 was obtained to 0{p'^), through a simultaneous fit of 
this parameter and the OPE corrections of dimensions six and eight to several 
spectral moments of the hadronic distribution. This determination is in good 
agreement with our 0{p^) result (120!) . Our quoted uncertainty has an smaller 
experimental contribution and includes a better assessment of the theoretical 
uncertainties. The value L^q = (—5.8 ± 0.2) • 10"^ (3.2 a smaller than ours) 
was extracted from r data in ref. [19] using the first "duality point" of the 
WSRs. The difference comes from underestimated theoretical uncertainties 
in this reference, as can be easily seen by choosing instead the second duality 
point or varying slightly the value of the first duality point. In fact the same 
reference [19] (see Eq. (10) therein) presents also a different estimate of L^q 
that is in very good agreement with our result. In ref. [40] both L^q and 
were determined, in good agreement with our findings which use the most 
recent 2005 data. An updated value of L^q, using the 2005 data, has also 
been given in ref. [21]. 

Our determinations of L\q{ij) and Cgj{n) at = Mp agree within errors 
with the \a.Yge-Nc estimates based on lowest-meson dominance [31,36,41,47]: 



- ~4Mj + 4M1 ^ "SMf ^ ' 



(23) 



(24) 



Eq. ( l22l) is also in good agreement with the result of ref. [38] for Csr based on 
Pade Approximants. These predictions, however, are unable to fix the scale 
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dependence which is of higher-order in I/Nq- More recently, the resonance 
chiral theory Lagrangian [36, 48] has been used to analyse the correlator 
n(s) at NLO order in the 1/Nc expansion [39]. Matching the effective field 
theory description with the short-distance QCD behaviour, the two LECs 
are determined, keeping full control of their yU dependence. The theoretically 
predicted values L'[o(Mp) = -(4.4±0.9)-10-=^ and Q7(Mp) = (3.6±1.3)-10-3 
GeV~^ [39] are in perfect agreement with our determinations, although less 
precise. A recent lattice estimate [49] finds L^q(Mp) = —(5.2 ± 0.5) ■ 10~^ at 
which is also in good agreement with our 0{p'^) result in (!20|) . 
A recent reanalysis of the decay vr^ e^vy [44], using new experimental 
data, has provided quite accurate values for the combination of Oij)^) LECs 
Lg + Lio. To one finds Ll{Mp) + L^^iMp) = (1.32 ± 0.14) ■ 10-^ while 

the 0(/) result Ll{Mp) + L[q{Mp) = (1.44 ± 0.08) ■ 10"^ is slightly more 
precise [44]. Combining these numbers with our results for L^q(Mp), one 
obtains 

f (6.54 ± 0.15) ■ 10-3, O(p^), 
LUMp) = I (25) 
[ (5.50 ±0.40) -10-3, C(/), 

in perfect agreement with the O(p^) result Ll{Mp) = (6.9 ± 0.7) ■ 10^^ of 
ref. [29] and the C(/) result Ll{Mp) = (5.93 ± 0.43) ■ lO'^ of ref. [45]. This 
last comparison represents an indirect check (in fact the only possible one 
for the moment) of our 0{p^) result for Liq. 

We have also determined the corresponding LEC of Liq in the SU(2) 
effective theory, both at LO and NLO: 

f 13.30 ±0.11, 0(p^), 
h= I (26) 
( 12.24 ±0.21, C(/). 

/^From a phenomenological analysis of the radiative decay vr Iw-f within 
SU(2) xPT, the authors of ref. [50] obtained ~Iq -'h = 2.57 ± 0.35 at C(p^), 
and iQ — h = 2.98±0.33 at 0{p^). Using these results and our determinations 
for I5 in ( l26l) . one gets 

15.87 ±0.37, C(p^), 

(27) 

15.22 ±0.39, 0{p^). 

At 0{p^) the comparison of these estimates of SU(2) LECs with previous 
results is straightforward, since they are proportional to the corresponding 
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SU(3) couplings, that we have aheady discussed. Our determination of I5 is 
the first one obtained at 0{p^), whereas for Ig ref. [51] finds Is = 16.0 ± 0.5 ± 
0.7, where the last error is purely theoretical, in good agreement with ours, 
although less precise. 
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